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Abstract 

We address the issue of the distribution of firm size. To this end we 
propose a model of firms in a closed, conserved economy populated 
with zero-intelligence agents who continuously move from one firm to 
another. We then analyze the size distribution and related statistics 
obtained from the model. Our ultimate goal is to reproduce the well 
known statistical features obtained from the panel study of the firms 
i.e., the power law in size (in terms of income and/or employment), 
the Laplace distribution in the growth rates and the slowly declining 
standard deviation of the growth rates conditional on the firm size. 
First, we show that the model generalizes the usual kinetic exchange 
models with binary interaction to interactions between an arbitrary 
number of agents. When the number of interacting agents is in the 
order of the system itself, it is possible to decouple the model. We 
provide some exact results on the distributions. Our model easily 
reproduces the power law. The fluctuations in the growth rate falls 
with increasing size following a power law (with an exponent 1 whereas 
the data suggests that the exponent is around 1/6). However, the 
distribution of the difference of the firm-size in this model has Laplace 
distribution whereas the real data suggests that the difference of the 
log sizes has the same distribution. 
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1 Introduction 



It is long known that the size distribution of the firms has a long tail which 
is remarkably robust [1]. A few very large firms can operate side by side 
with a large number of small firms. Ref. [2] presents clear evidence that the 
distribution can be characterized very well by a power law and regarding the 
stability of the law, as the same reference puts it, this feature has survived 
changes in the political, regulatory and social regimes; the last one being 
caused by the demographic changes in the work force due to the influx of 
women in the labor force. Also, numerous innovations and technological 
changes in the production process were unable to affect it. Lastly, firm 
mergers, acquisitions, death and birth of firms did not affect this feature. 
Ref. [3] studies the panel data on the firm dynamics (all publicly traded 
US manufacturing firms in the time span 1975-1991) and concludes that the 
growth rates show two more significant features. One, the distribution of the 
growth rates of the firms has an exponential form and the standard deviation 
of the growth rates of the firms fall with increasing firm sizes following yet 
another power law. 

The above finding indicates that the statistical features for the firm growth 
process are independent of microeconomic decision-making processes (at least, 
to a first approximation) like why people choose to leave their job etc. Hence, 
we do not indulge in providing any microeconomic foundation for the firm 
dynamics. However, the rate at which the firms gain and lose workers is of 
interest to us. This rate is called the turnover rate in the economics litera- 
ture. An alternative but closely related interpretation of the turnover rate is 
that it measures how long the employees stay in their respective jobs. It may 
be noted that the rates of hiring and separation for developed economies are 
very high. For example, in USA (2009-2011), the average total seasonally 
adjusted annual hiring rate and separation rate was around 38-40% (see Ref. 
[1]). Hence, the turnover rates (interpreted as the average length of employ- 
ment) may be very low. We intend to show that the turnover rates play a 
crucial role in the firm size distribution and related issues. One important 
aspect of job separation and worker hiring is that the process follows the 
rule of local conservation. If one worker goes from one firm to another then 
the total workforce remains unchanged but the workers' distribution across 
the firms change. Since the workers at any given year (or quarter) move 
around in a very large number of firms, we model this process as a repeated 
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interaction between a large number of agents (firms) which exchanges a finite 
amount of (number of) workers between themselves. Clearly, the idea of the 
kinetic exchange model is suitable for this purpose. 

In this paper, we present a model of firms in a closed, conserved economy 
populated with zero-intelligence agents. The firms are modeled as collections 
of agents who continuously move from one firm to another. But on an aggre- 
gate level, there is no fluctuation. The firm's size is solely determined by the 
number of agents working in the firm. We also assume that time is discrete. 
The basic idea of the model is that each period there is a group of agents 
in each firm who wants to move to another firm with some personal motives 
(like utility maximization due to wage increase etc.) which we do not model 
here. We call the turnover rate A. Hence, (1 — A) fraction of each firm's 
workforce would want to move out of the respective firms. Each period there 
will be a pool of such agents who wants to shift from one firm to another. 
Some of them will find a new job (that is they will move to new firms) while 
the rest has to continue in their earlier position. There is no unemployment 
in the model. This process is repeated until the distribution of workers settle 
down and this distribution will be the size distribution of the firms. We shall 
show that this model is closely related to the kinetic exchange models of 
markets and in fact, it generalizes the usual binary trading (collision) model 
to interactions between an arbitrary number of agents (firms in this case). 
The exact distributions in some cases, will be provided. Subsequently, we 
shall study a modification of the basic process which leads to a Power law in 
the size distribution of the firms. Then we show that the standard deviation 
decreases with increasing size following yet another power law and we study 
the corresponding distributions of the growth rates. 

The related literature is varied and vast. Ref. [1] is probably the first sys- 
tematic treatment of the subject. To model the growth of a firm it suggested 
a stochastic process which essentially states that the growth rate of the firm 
is independent of its size. This prediction and its result that the distribution 
of the firm size would be log-normal, was found to be approximately correct 
[5]. However, there are evidences that the formulation was not entirely cor- 
rect. Ref. [H] first observed that the standard deviation in the growth rate 
falls as the firm size increases. This finding is supported in later studies as 
well (see Ref. [31 [TJ El [U [10]). Ref. [11] collects data for US manufacturing 
firms (1974-1993) which support their earlier finding in The more im- 
portant part of this study is that they showed that the proposed statistical 
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features are robust to firm birtli and death due to merger or bankruptcy. For 
an overview, see Ref. [12] and references therein. See also Ref. [13] for a 
detailed analysis of the Gibrat's law, Pareto index and Pareto law. The data 
set used is mainly the panel data of the Japanese firms. On the theoretical 
ground, Ref. [H] contains study on stochastic properties of the dynamics 
of firm growth. Ref. [12] presents a model of hierarchical organizations to 
explain the observed regularities. See Ref. [Ml |T7] for a separate theoretical 
approaches to the dynamics of company growth. See Ref. [TH [191 1201 [21] for 
detailed discussion of the kinetic exchange models of markets. Lastly, Ref . 
[22] provides theories of job matching and turnover. 

This paper is organized as follows. In section [21 we propose the basic 
model and derive the exact distribution of workers in the firms. In the next 
section, we modify the model which leads to the power law distribution of 
the firm sizes. In section [H we show the distributions of the growth rates of 
the firms in this economy and compare them with real data. Then follows a 
summary and in the appendix, we have presented short discussions on how 
the kinetic exchange models are related to the generalized Lotka-Volterra 
equations and also, how can we apply the model stated in this paper directly 
to model the income/wealth distributions. . 

2 The model with constant turnover rate, A 

We assume that time is discrete. The economy consists of an array of N firms 
which can absorb any amount of workers that come to it. The workers are 
treated as a continuous variable (infinitely divisible). At the very beginning 
of the process, all firms have exactly one unit of workers (more formally, the 
measure of workers is one for each firm). The fraction of workers that decides 
to stay back in their firm (which we interpret as the turnover rate), is denoted 
by A which may vary between the firms. For the time being, we treat them as 
given and constant across the firms. This treatment is pioneered by Ref. [23] 
in the context of modeling income/wealth distributions. As we said earlier, 
the firm's size is just the measure of workers working in the firm. There are 
other indicators of the firm size as well e.g. quantity of goods produced, sales, 
cost of goods sold, assets or value of the properties. But we note that not all 
firms produce the same, identical goods. The inputs also differ very much. 
To consider capital holding (or the value of assets), that is not always easy 
to measure (large fluctuations happen in the stock market in short spans of 
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time adding or wiping huge amounts from the value). Hence, it is easier to 
use the size of the workforce as a proxy for the firm size. We denote the 
firm size of the i-th firm (that is the work force) by Wi {i < N where N 
is the number of firms). Also, suppose that the number of firms from which 
the workers are leaving and moving into, is n. At each time point (1 — A) 
fraction of the workforce of those n firms wants to leave. So there would be 
a total pool of workers that wants to change their workplace. Next, this pool 
of workers is randomly divided into those n firms. Hence, the dynamics is 
given by the following set of equations. 



wi{t + l) = A^;i(t) + ei(t+i)(l-A)^it;,(t) 



Wi{t + 1) = Awi(t) + ei(t+i)(l - A)^w;j(t) 



Wnit + 1) = \Wn{t) + e„(t+i)(l - X)'^Wj{t) (1) 

such that ^j(t) = 1 for all t. As is evident from above, this is a straight 
generalization of the usual kinetic exchange models (with n = 2) that has 
primarily been used to study the income/wealth distribution models (see Ref. 
[TSl [201 EI] ) • Regarding the notations, we use t within the first bracket when 
referring to the endogenous variables like the size of the firm ( w{t)) and we 
use the same in subscript when referring to the exogenous random variables 
(e.g., e^). Similarly, with a slight abuse of notation, we denote the probability 
density function (pdf) of the exogenous random variable x by f{x). However, 
for the distributions of the endogenous random variables (for example, the 
distribution of the firm size or of the growth-rate), we use P{.) {P{w) and 
P{g) resp.). 

Construction of e 

Here, we consider a few constraints on e. 

1. The sum of all e^s has to be equal to one. 

2. The expectation, E{ei) = 1/n for all i and the distributions of all are 
identical. 
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3. If n = 2, ej ~ uniform[0, 1]. We impose this constraint so that at the 
lower hmit of n, we get back the usual CC-CCM models (see Ref. [TS]). 

Formally, the problem then boils down to that of sampling uniformly from 
the unite simplex (see Ref. [21] )• We follow the standard algorithm and 
below we derive the distribution of e. 

1. Create a vector of independent random variables drawn from uniform 
distribution over [0, 1], ^i, ^2, • • • , ^n- 

2. Take logarithm of all the elements of the vector and multiply the ele- 
ments by -1. 

3. Divide each element by the sum of all the elements. Call the i-th result 
e, for all i. 

We derive the probability density function of the below. Consider ei first 
for simplicity. The probability that ei is less than some 6 is 



Prob.{ei<9) = Pro6.(- In ^1 < -^ln(ei6---en)) 

= l-Pro6.(ei<#e|^...e#') 
/•I _0_ 

= 1 ^ [ / ^2 "'^^2]""^ (using independence) 



^0 

= l-(l-^)"-\ 

This is true for all e^. Therefore, the pdf of is 

/(6,) = (n-l)(l-e,)"-^ (2) 

that is, e has a beta pdf with parameters 1 and n — 1. Clearly, when n = 2 
the distribution of e is uniform [0, 1] as expected. 

2.1 Reduced form of the model 

First, we note that the solution to the usual kinetic exchange model with 
binary interaction is not known yet (see Ref. fiSi [19] ). The resultant distri- 
bution is approximated by gamma probability distribution function [25] . But 
moment considerations show that the distribution does not have a gamma 
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Figure 1: 

Workers' distribution across the firms: comparisons between binary 
interaction i.e., usual kinetic exchange model (blue squares; see Ref. [T8] ) 
and A^-ary interaction model (red stars). Three cases are shown above, viz., 
A = 1/4, A = 2/4, A = 3/4. All simulations are done for O(IO^) time steps 
with 1000 agents and averaged over O(IO^) time steps. Note that in Sec. 
12.21 we have derived that for A = 0, both curves are identical. Discrepancies 
appear for A > as is apparent in the above diagram. 
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form (see Ref. [20]). Here, we derive an exact result of the case where the 
number of interacting firms is in the order of the system size N i.e., we 
consider the case where 2 « n < N . 

Note that if n is of the order of A^, Wj is well approximated by n (recall 

that E{wj) = 1 for all j). To make sure, note that ^ by the 

structure of the model. For exactness, we shall assume that all firms interact 
at every step, i.e., n = N. Then the system of equation becomes 



wi{t + l) = Awi(t) + ei(t+i)(l-A)iV 
Wi{t + 1) = Xwi{t) + ei^t+i){l - X)N 



Wnit + 1) = XWnit) + enit+l)il - X)N (3) 

with each having a beta distribution as has been found in Eqn. [2] (see 
Construction of e in Sec. [2]). Note that in this form, we get rid of the 
effects of Wj{t) in the evolution equation of Wi{t) for all j ^ i. One more 
simplification is possible. Let fi = N{1 — A)e ignoring the subscripts. Given 
A^, it is easy to verify that the probability distribution of /i is 

Hence, for large we can approximate the distribution as the following, 

^im /(/i) ~ ^e~'^'^ where ip = (5) 

Therefore, the system reduces to 

wi{t + i) = Xwi{t) + fim+i) 



Wi{t + 1) = XWi{t) + /i2(t+l) 



WN{t + l) = XwN{t) + fJ^Nit+l), (6) 

which is a system of autoregressive type equations with the distribution of 
errors (/i) given by Eqn. [5l 
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2.2 Steady state distributions 



2.2.1 For A = 

In the above section (Sec. I2.ip . we have derived the reduced form equations. 
Below, we find their solutions. One noteworthy feature is that with A = 0, 
the system further reduces to 

Wi(t) = Neu for all i. (7) 

As we showed in the above section, the steady state distribution would be 
exponential. Note that the result is identical to the case where the interaction 
is binary. Also, we can provide another proof by conjecture. Let us rewrite 
the system as 

Writ) = eu Wj{t - 1) j for all z. (8) 

Let us conjecture that the steady state distribution is exponential. More 
precisely, let f{wj) = exp{—Wj) for all j. Clearly, J2f "^ji^ — 1) has a gamma 
pdf with parameters 1 and A^. Recall that e has a beta pdf with parameters 
1 and A^ — 1. Therefore the distribution of their product is again exponential 
(see Thm. 2.3 in Ref. [27]) confirming our conjecture. 



2.2.2 With positive A 

However, the above result (the equivalence between the distributions gener- 
ated by binary interactions and N-ary interactions) does not hold in presence 
of positive A. First, we discuss the discrepancies in the second moment. Then 
we move on to derive the exact distribution. We denote the central moment 
of order n > 1 of a variable 

E{x - E{xY) = E{J2 ( ^ ) x'E(-x)("-')). 

For n = 2, E{x — E{x)"') corresponds to the variance of x and is denoted by 
V{x). Since the system is conservative and the initial workforces (i.e., the 
firm sizes) were unity for all firms, it is obvious that E{wi) would be unity. 
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So we can write the n-th moment of the distribution of size without subscript 

as 

E((^-in = i?(X:(^)(-^)')- (9) 

We assume that Wi and Wj are independent variables (technically, they are 
not since the sum of all WiS is constant, in this case; but for large this 
is a good approximation ). It is easy to verify that with all firms interacting 
(n = A^), the variance is given by 

whereas in the case of binary interaction [21] 

(1-A) 



V{w) 



(1 + 2A)' 



Note that for A = 0, variance is unity in both cases which is consistent with 
our derivation that the distribution is the same (exponential) in both cases. 
Let us write the system as 

w{t + 1) = Xw{t) + fit+i 

which can be rewritten with the lag operator L as (1 — XL)w(t) = fit and 
hence, 

w{t) = /it + A/it_i + AVt-2 + AVt-3 + • • ■ • (10) 
Recall that (Eqn. E> 

Therefore in the steady state, 

W = flo + fl2 + fi3 + ■ ■ ■ ■ (11) 
where flj is distributed as 

1 
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We can neglect the terms with high powers (more than say k) of A. Then w 
is essentially the sum of k exponentially distributed random variables with 
different parameters. Note that the Laplace transformation L{s) of Hj is 
4>j/{.(f>j + s) with = 1/(A-^(1 — A)). Since the fi/s are i.i.d., pdf of w would 
be the convolution of the pdfs of the k random variables. By property of 
Laplace transformation, it can be verified that the distribution of w would 
be (see Ref. |2B] for detailed discussions and different proofs) 

f{w) = Y,<l>^exp{-<Piw) n (t^) (12) 

with 0j defined as 0j = 1/(A*(1 — A)) (see figured]). 



3 Distributed turnover rates, ^ Xj 

So far, we have considered only fixed A. In this section we consider dis- 
tributed A (i.e., the turnover rates differ across firms but they are fixed over 
time) following Ref. |T8]. Specifically, we assume that the turnover rates 
are uniformly distributed over the interval [0, 1] across the firms. The new 
system of equation is 

n 

Wi{t + 1) = Aii(7i(t) + ei(t+i)^(l - Aj)wj(t) 

j 



Wi{t + 1) = \iWi{t) + ej(t+i)^(l - \j)wj(t) 

j 



Wn{t + 1) = A„w„(t) + e„(t+i)^(l - Xj)wj{t) (13) 

i 

To solve Eqn. [iniin the steady state, note that {l — \i)E{wi) = C, a constant, 
solves the problem. Therefore, we can rewrite the system of equation as 
(assuming n = N) 
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Figure 2: 

Finding the value of C from Eqn. [TH We have considered three system 

sizes viz. = 100 (uppermost), 200 (middle) and 300 (lowermost). 
Clearly, the value of the constant C decreases with increasing system size 

N . See also figure [31 
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Figure 3: 

Upper panel: Dependence of C on the system size A^. As we can see the 
value of C falls rapidly with increasing system size (for small systems). 
With N = 3000, C ~ 0.1285. Lower panel: Dependence of C on l/log{N). 
We simulated systems with different sizes 
(N = 5000, 10000, 15000, . . . , 60, 000) for - 10^ time periods. The 
observation fits well with 0.1137 + 0.095/log{N). The rightmost point 
(iV = 5000) is above the fitted line because of the effect of small system 

size. 
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Wi{t+1) 



Wiit + 1) = \iWi{t) + C/i2(t+l) 



WNit+l) = XNWN{t) + C^Nit+1), (14) 

with Aj ~ uniform[0, 1]. Also recall that /i, = iVe, with /(/ii) = exp{—^i) 
(see Sec. I2.2.1I) . Ref. ^29j finds the value of the constant C, in the context 
of the usual kinetic exchange models with binary trading scheme. Here, we 
confirm by simulation that {1 — Xi)E{wi) is actually a constant for all i (given 
the system size i.e., N; see figures [2] and [3]) . Hence, we can regard Eqn. [H] 
as correctly representing the model (see also Ref. |30j which treated markets 
populated with agents each having a different autoregressive process defining 
their wealth evolution). The resultant distribution of the above model is a 
power law (see figure Hj). Ref. [31j considered a slightly more general version 
of this type of maps. Following Ref. [29] , a very simple proof is considered 
below. Note that (in the steady state) by taking expectations on both sides 
of Eqn. [131 we can rewrite it as 

(1 - Xi)E{wi) = C. (15) 

By taking total differentiation and rearranging terms, we get 

dX _2 
dw ' 

where w represents E{w). Hence, the average workforce in a firm with a 
particular A is given by Eqn. [151 Also, the relation between the distribution 
of A (i.e., /(A)) with that of w is given by the following Eqn. 

P{w)dw = f{X)dX. 

The last two equations show that in an array of firms with uniformly dis- 
tributed A, the distribution of w would be 

Hence, the firm size has a power law distribution (Zipf's law; see Ref. [2]). 
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4 Growth rate 



Truly speaking, in this model there is no absolute growth. The economy 
as a whole is completely conserved. There are a few firms (with high A) 
that grows initially. But after achieving their average size, they do not grow 
any further in the absolute sense. But of course, fluctuation is still present. 
Recall that the reduced form of the model is given by 

Wi{t + 1) = XiWi{t) + Ciiiit+i) for all z, (16) 

which is simply a set of autoregressive type equations. Let us define growth 
rate as rj = Wiit + 1) / Wiit) . Clearly, rj = Aj + C fii(t+i) / Wi(t) in this model. 
Evidently, as A, rises, the average size E{wi) also rises. Therefore, the vari- 
ance (or the standard deviation) of the growth rate falls with increasing size 
(see figure [5]). This model captures this feature well though it does not match 
the exact exponent. The standard deviation decreases following a power law 
with exponent -1 (see figure [H]) whereas the real data set suggests that the 
exponent is 0.16 ± 0.03 (see Ref. [3]). We also studied the distribution of 
the growth rate of the firms. Ref. [3] defined growth-rate as logr with r 
defined as above. However, we can approximate the growth rate as follow- 
ing logrit = log{wi{t) / Wi{t - 1)) ~ {wi{t) - Wi{t - l))/wi{t) (by adding and 
subtracting 1 to r). Note that for A = 0, the distribution of any firm (i.e., 
w{t)) would be exponential. Therefore, the numerator in the expression of 
logrj has a Laplace distribution (see figure [7]). But the growth-rate logr as 
has been defined in [3] clearly does not have a Laplace distribution in this 
model. In fact, in no way it resembles the proposed Laplace distribution (it 
has too many discrete jumps, specially for firms with small A ). Hence, our 
model does not perform well to reproduce the pdf of growth rate found in 
Ref. [g. 

5 Summary 

We have studied a model of firm dynamics. There are a number of well known 
statistical regularities in the firm dynamics (see for example Ref. [21 [3]). The 
main features considered here are (a) power law decay in size distribution, 
(b) reduction in fluctuation in growth rate with increasing size of the firm 
(following another power law) and (c) Laplace distribution of the growth rate. 
Ref. mentions another regularity concerning the distribution of payments 
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Figure 5: 

Standard deviation of rj(t+i) = Wi(t + l)/wi(t) decreases with increasing A 
(upper panel) and average size (lower panel). It shows exponential decay 
with respect to A for < A < 0.9. The straight line shown in the upper 
panel is 5.exp{—AX). All simulations are done for 0(10'') time steps with 
5000 agents and averaged over O(IO^) time steps. 
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Figure 6: 

Standard deviation of \ogri(^t+i) = ^og Wi(t + 1) — log Wi(t) decreases with 
increasing A (upper panel) and average size (lower panel). Clearly, it shows 
a power law decay with respect to size as has been documented in Ref. [3] . 
However, the exponent found in this model is -1 whereas data suggests that 
it is 0.16 ± 0.03. All simulations are done for O(IO^) time steps with 5000 
agents and averaged over O(IO^) time steps. 
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Figure 7: 

Distribution of growth rates gi{t) = Wi{t) —Wi{t — 1) for turnover rates Xi = 
0, 0.2, 0.4, 0.6, 0.8, 0.95. Evidently for small Xi, gt has Laplace distribution 

(bi-exponential) . As lim A — )■ 1, the distribution becomes one sided 
exponential only (see Eqn. [16] and note that the error term is exponentially 
distributed). However, contrary to what we get here, the data suggests that 

\ogWi{t) — \ogWiit — 1) has a Laplace distribution (Ref. [3J). All 
simulations are done for O(IO^) time steps with 5000 agents and averaged 

over O(IO^) time steps. 
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to the workers in different firms. But we liave neglected it completely because 
we did not consider any strategic behavior on the part of the workers or the 
firms (i.e., the firm owners) . 

The model that we propose is a multi-agent model with n-ary interactions 
{2 < n < N) at each time step. We present some analytical results on the 
steady state distributions for constant turnover rate (across the firms) ob- 
tained from the model which seems to not agree with the earlier approximate 
results (see Ref. [25]). Specifically, we show that if the number of interacting 
firms goes to infinity the exact distribution of the firm size is given by an in- 
finite sum of weighted exponentials. Hence, at least in this limit the solution 
is definitely not gamma distribution as had been proposed for binary interac- 
tion. With distributed turnover rates, the model can very easily produce the 
power law in size distribution as is done by the usual kinetic exchange models 
with binary interactions. However, the importance of the generalization from 
binary to n-ary interaction lies in the idea that it better captures the work- 
ers' flow among a huge number of firms. Next, we study the distributions 
of the growth rates of the firms. We show that this model quantitatively 
captures the observation that the fluctuations in growth rates fall according 
to a power law with increasing firm sizes, though it fails to match the exact 
exponent (the model gives 1 whereas the data suggests 1/6). Another short- 
coming of the model is that, in this model, all growths are relative i.e., the 
economy as a whole is not growing which is certainly not the case with the 
real economies. The major point of departure of our model from The law of 
proportionate effect (see Ref. [IjJ which assumes zero effect of the firm size 
on its growth rate) is that we assume autocorrelation in the growth rate. In 
fact, we present the firm's growth process by an autoregressive process of 
order one so that the growth rate is affected by the size. See Ref. [22] which 
documents autocorrelations in the firm growth processes. Ref. [53] claims 
that small firms actually show a negative autocorrelation whereas the large 
firms have positive autocorrelation. However, the evidence is not conclusive. 
In the appendix, we discuss briefly about how this generalized version of 
the kinetic exchange models are related to the generalized Lotka-Volterra 
model and also, how can we apply the model stated above to model the 
income/wealth distributions. 

Acknowledgement: The author is grateful to Bikas K. Chakrabarti and 
Arnab Chatterjee For some useful comments. 
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6 Appendix 



6.1 Comparison with the generaUzed Lotka-Volterra 
(GLV) model 

Ref. [35j raised a question that whether and how, the kinetic exchange models 
are related to the generalized Lotka-Volterra model (both can produce power 
law distributions). One major obstacle in answering the question was the 
problem that the kinetic exchange models focused on binary interaction and 
GLV takes in to account all the agents in each interaction. However, given 
the above formulation, one can directly compare the two mechanisms. 

Ref. [31] presents the GLV mechanism by the following equation (for 
l<i<N), 

Wiit + 1) = \{t + l)wi{t) + a{t)w{t) - c{t)wi{t)w{t). (17) 

Two notable differences of GLV with the model proposed above are the pres- 
ence of a time varying A and a nonlinear interaction via the average w {w) in 
the GLV. Recall that in the model proposed as a generalized kinetic exchange 
model, the average is always fixed (at unity, in this case). Ref. [31] reduces 
the system to N decoupled equation of the following form, 

ik{t + l) = X{t)vi{t)+a{t). (18) 

Comparing Eqn. [18] with Eqn. [161 we see that the essential difference be- 
tween the two systems is whether A varies over time or not. Another impor- 
tant point is that in GLV, A has to be greater than 1 sometimes which is not 
possible in the other case. In short, we can say that the GLV mechanism de- 
pends on the process of random multiplicative maps whereas the generalized 
kinetic exchange model does not. 

6.2 Generahzed exchange model as an equilibrium out- 
come 

In this subsection, we discuss how the model with A^-ary interaction be ap- 
plied to model wealth/income distribution. Though we think that binary 
trading is much more common in the real market place (than A^-ary trad- 
ing), we present briefly a direct generalization of the framework presented in 
which dealt with binary trading mechanism. 
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Let there be agents each having 1 unit of perfectly divisible money in 
their possession (at the beginning of all trading). At each period each of 
them produces Q unit of commodities {Qi may be different from Qj for all 
i and j) such that no two commodities are the same. Let the preference of 
the i-th agent be defined as 

TJ. — ^"1^02 ^OLf^ am 

LJ I — Uy-^ ^2 . . . -^TY ' ' ''i 

The budget constraint would be piXi + P2X2 + . . . + nii < Mi + PiQi. We 
make the standard assumption that ai + 02 + ■ ■ ■ + X = 1 where A = 
is the savings propensity. Then we can write the constrained maximization 
problem as 

C = x1^X2^ . . . - Ji{piXi + P2X2 + . . .nii - Mi+ PiQi) 

where p, is the Lagrange multiplier. Solving the optimality conditions, we 
get (denoting by A) 

m* = \{Mi+piQ{) 
for all i. Solving for equilibrium price vector, one derives 

mi{t + 1) = \mi{t) + 6^(1 - A) y^mj(t) 

where can be suitably defined as a beta r.v.. Clearly, the above equation 
reduces to the following 

mi{t + 1) = \mi{t) + 
(this is exactly the system we studied above in section [2]). 
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